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METHOD OF COMBINATORIAL MULTIMODAL OPTIMISATION 

The present invention relates to methods of combinatorial multimodal 
optimisation, and in particular to the use of genetic algorithms for solving 
5 combinatorial multimodal optimisation problems. The invention is expected to 
find particular although not exclusive application in the fields of load balancing 
(e.g processor and network load balancing), scheduling optimization (including 
production scheduling, resource assignments, timetable scheduling) and 
resource planning. 



An optimisation problem is always a simplication of a real world problem, an 
engineer needs to suggest various possible optimal solution alternatives to the 
decision maker in order to understand this problem better. 

15 Many problems in the scheduling of manufacturing systems and artificial 
intelligence can be regarded as combinatorial multimodal optimisation 
problems, which require the determination of multiple integer solution vectors 
that maximise (or minimise) a given objective function with regard to some 
definite set of constraints. Problems of this class, which are characterised by 

20 their discreteness, high dimension and noise, are quite different from 
multimodal continuous function optimisation problems. The combinatorial 
multimodal optimisation problem is arguably one of the most deceptive 
problems for GA applications. 



10 



25 So-called genetic algorithms (GA) are well-known stochastic optimisation 
search methods which are based on a global search procedure. But the use 
simple genetic operators (selection, multi-point crossover and mutation) to 
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solve a typical combinatorial multimodal optimisation problem, provides only a 
single (global or local) optimal solution for each search procedure [He Liwen 
and Mort. Neil: Genetic Algorithm for Scheduling Optimization on a Multi- 
Product Batch Processing Machine. 4 th IFAC Workshop on Intelligent 

5 Manufacturing Systems, Seoul, Korea, PP165-170, 1997]. The reason is that, as 
there is no selective pressure among any of the peaks for the finite population in 
the fitness landscape, the population converges to one alternative or another 
randomly. This problem is known as genetic drift ~ stochastic errors in 
sampling caused by small population size. In other words, this difficulty comes 

10 from the fundamental theorem of genetic algorithms which tells us that 
exponentially increasing trials will be given to the best schemata observed, thus 
losing some potentially useful schemata during the evolution. 

Therefore, the problem is how to reduce the effect of this genetic drift and keep 
15 the population diversity in GA search in order to find different fitness peaks, 
and lead populations to move from local optima to global optima. 

A typical NP-hard Combinatorial Optimisation Problem could be described as 
follows [Jelasity, M. and Dombi, J: GAS, a Concept on Modelling Species in 
20 Genetic Algorithms. Artificial Intelligence, Vol. 99, No. 1, pp. 1-19, 1998]: 

A set W= {w u w 2 , ... , w n } of n integers and a large integer P is given. We 
would like to find anScfF such that the sum of the elements in S is close to, 
without exceeding, P, perhaps subject to some problem-specific constraint 
conditions. The number of elements in S is given as m. 

25 

To take a simple example, let us assume that we have eleven jobs of various 
sizes that need to be shared, within a computer system, between five parallel 
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processors. Each processor is capable of handling a single job of size 9, or any 
combination of smaller jobs having a total aggregate size of more than 9. 

In this example, let us assume that the eleven jobs are defined by the set W = 
5 {2, 7, 5, 4, 3, 6, 2, 4, 2, 8, 1} of 11 integers; using the above terminology, n = 
1 1, m = 5, and P = 9. The task, then, is defined 5 groups of jobs such that the 
sum of the values within each group is close to, but does not exceed, 9. 

One possible solution is as follows: 

Processor 1: {2,7} 
Processor 2: {5, 4} 
Processor 3: {3, 6} 
Processor 4: {2,4,2} and 
Processor 5: {8, 1} 

With this solution, four of the processors have an aggregate workload of 9, 
while the remaining processor (processor 4) has an aggregate work value of 8. 
Since all of these workloads are sufficiently close to but not exceed the 
20 maximum value of 9, we say that the fitness value of this possible solution is 5. 
Since that is the same as the number of processors, this solution is optimal. 

Other trial solutions may give less good results. In the following groupings, 
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Processor 1: {2,7} 
Processor 2: {5,4} 
Processor 3: {3,6,2} 
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Processor 4: {4,2}; and 
Processor 5: {8, 1} 

only the first, second and fifth processors have a full workload. The third 
5 processor has a workload of 11 and is so overextended, while the fourth 
processor has a workload of 6 and is under-utilised. Since only three of the 
five groupings meet the problem criteria, the fitness value of this possible 
solution is 3. 

10 If the number of the global optima is greater than 1 in the solution space, then 
this problem is regarded as combinatorial multimodal optimisation problem. 

Conventional genetic algorithms start with an initial population of individuals 
(trial solutions). The fitness value of each individual is evaluated according to 

15 some function (sometimes called the "objective function"), and there then 
follows a solution process under which individuals with higher fitness values 
have a higher probability of surviving into the next generation. Typically, 
individuals with low selection values are eliminated entirely, and individuals 
with high selection values are replicated several times, in proportion to their 

20 fitness value. Genetic operators such as crossover and mutation are then 
applied to the population, and the fitness values of the new individuals are re- 
evaluated. This process continues until the algorithm converges. 

As mentioned previously, such an approach is limited because of the problems 
25 of genetic drift, and the fact that only a single optimal solution can be obtained. 
Indeed, the very nature of the procedure tends to ensure that sub-optimal 
individuals (in terms of fitness value) ultimately become reduced and 
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eliminated from the population, so ensuring that as the algorithm proceeds only 
a smaller and smaller area of the solution space is explored. That may be 
useful for finding local optima, but at the expense of finding global optima. 

5 Several approached have been proposed to deal with these problems, such as 
niching methods, fitness sharing, [Sareni, B and Krahenbuhl,. L: Fitness 
Sharing and Niching Methods Revisited. IEEE Transactions on Evolutionary 
Computation. Vol. 2, No. 3, 1998] and parallel genetic algorithms [ 
Chipperfield A.J and Fleming, P.J: Parallel Genetic Algorithms: A Survey. 

10 Research Report No. 518, Department of Automatic Control and Systems 
Engineering, University of Sheffield, 1994]. Other techniques include those 
proposed by Sirinivas, M. and Patnaik. L.M: Adaptive Probabilities of 
Crossover and Mutation in Genetic Algirithms. IEEE Transaction on Systems, 
Man and Cybernetics. Vol. 24, No. 4, 1994; Schneider. G., Schuchhardt. J. and 

15 Wrede. P: Evolutionary Optimization in Multimodal Search Space. 
Biological Cybernetics. Vol. 74, No. 3, pp.203-207, 1996; Yao. X and Liu. Y: 
Fast evolution strategies. Control and Cybernetics. Vol.26, No. 3, pp.467-496, 
1997; and Jelasity, M. and Dombi, J: GAS, a Concept on Modelling Species in 
Genetic Algorithms. Artificial Intelligence, Vol. 99, No. l,pp. 1-19, 1998. 

20 

Both niche techniques and adaptive GA require a previous understanding of the 
number and the distribution of fitness peaks, and need careful setting of various 
parameters. Parallel GA has proven very reliable and efficient in exploring 
multimodal fitness peaks, and it also needs fewer assumptions than niche 
25 methods. The main limitation of implementing Parallel GA is the sub- 
population size which limits the total number of global optima found in the 
search. 
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According to the present invention there is provided a method of combinatorial 
multimodal optimisation for finding multiple optimal ways of dividing a set W 
of n values into m groups, such that each of the groups satisfies a respective 
5 constraint condition, the method comprising: 

(a) defining an initial population of individuals, each representative 
of atrial solution; 

(b) calculating for each individual a fitness vector indicative of 
whether the constraint condition for each group has been satisfied; 

(c) selecting a plurality of individuals for the next generation in 
dependence upon their respective fitness vectors; 

(d) creating a new population including the selected individuals; and 

(e) repeating steps (b) to (d) until the population stabilizes, the 
individuals of the stable population representing multiple optional 
ways of dividing the set W. 

Experimental results show that both convergence speed and quality are greatly 
improved in comparison with other GA techniques. They also suggest that this 
20 novel method has the capability of generating sufficient genetic innovations, 
preserving promising genetic diversity, guiding population escape from local 
optima and exploring multiple global optima effectively and efficiently in this 
difficult combinatorial multimodal optimization problem. 
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The novel method works through deeply preserving the population diversity 
based on a phenotypic analysis, rather than simply a genotypic analysis. This 
reinforces the power of the method to process many building blocks in the 
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search space simultaneously. This power is stronger than that of the implicit 
parallelism found in simple genetic algorithms, and thus enhances the 
exploration of GA space for multiple optima. 

5 In this novel method a hybrid drift model, through the complementary 
interaction between genetic drift and neutral drift, can help a population escape 
away from local optima. When the relative fitness vector of one individual is 
non-dominated by that of any other one in the entire population, 'neutral drift' 
is in effect. Any non-dominated individuals are regarded as having equal- 

10 (relative)-fitness values and kept for the next generation, thus facilitating the 
exploration of the different adaptive peaks. On the other hand, when one 
individual is dominated by another, 'genetic drift' effectively lets this less fit 
individual die out. This model has been shown to be effective in leading 
populations to escape from local optima to global optima in this typical 

1 5 combinatorial multimodal optimisation problem. 

The invention may be carried into practice in a number of ways and one 
specific example will now be described, merely by way of example. 

20 Briefly stated, the preferred method of the present invention makes use of a 
genetic algorithm having the following steps: 

1. An initial population of individuals is randomly generated, each 
individual being represented by a string of integers which defines a trial 

25 solution; 

2. A fitness value and fitness vector is calculated for each individual, 
according to an objective function; 
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3. On the basis of the fitness values and the fitness vectors, suitable 
individuals are selected to form a sub-population whose size is a portion 
of the overall population; 

4. A roulette wheel selection method is men used to create offspring, from 
the selected individuals, to form a new population; 

5. Highly uniform crossover and mutation is applied to the individuals of 
the new population; and 

6. Steps 2 to 5 are repeated until the algorithm converges. 

Each of these individual steps will be described in more detail below. 

First, we will consider representation and initialisation of the population. A 
major difficulty in applying GA in a scheduling problem is finding an 
appropriate representation for the population. Here, an integer string 
representation of the individual is selected according to the following 
conditions: 

a. Completeness: the string representation should contain all possible 
individuals in the search population. 

b. Uniqueness: the string maps to the individual on a one to one basis. 

A second desideratum is to generate the initial population randomly. For 
example, a population is composed of N chrom individuals each of which is an 
integer string with n genes. The initial population is produced by creating a 
(Nchrom x ") matrix with the component randomly selected from the set {1, 
m}. 

In contrast with the prior art, the fitness of each individual is determined not 
only on the basis of the fitness value but on that individual's fitness vector. 
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The fitness vector is a binary vector of length m, the individual bits of which 
indicate, for that individual, whether the groupings satisfy or do not satisfy the 
problem constraints for each of the m groupings. 

5 To go back to the example previously discussed, in which 1 1 parcels of work 
have to be divided between five processors, such that no processor takes more 
than 9 work units, the first solution given satisfied the constraints for all of the 
five processors. Its fitness value was therefore 5, and its fitness vector 
FV = [1 1 1 1 1 ]. The second, less-good solution, satisfied the constraints only 
10 for the first, second and fifth processors, and hence had a fitness value of 3. 
The fitness vector in that case becomes [110 0 1]. 

It will be understood of course that, with this approach, the fitness value is 
simply equivalent to the number of positive bits within the fitness vector. 

15 

We now introduce some additional terminology, to be used in the remainder of 
the description. GA applications sometimes refer to the individuals within the 
evolving population as genotypes. Following that analogy, we shall refer to 
analysis using the fitness vector as phenotype or phenotypic analysis. We 
20 define the groups of combinations making up the solutions, for each individual, 
as a phenotype schema. 

To clarify this further, it may be useful to consider a further example which 
comes from a practical shoe batch processing machine (BPM) optimization 
25 problem. This problem is described in Fanti, M.P, Maione, B, Piscitelli, G and 
Turchiano, B: Heuristic scheduling of jobs on a multi-product batch processing 
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machine, International Journal of Production Research, 34, pp. 2163-2186, 
1997. In this problem, n is 21 , m is 7, P is 8 1 1 , and the vector Wis defined by: 

[12 12 170 176 216231 266 310 380 497 714 144 128 153 660 688 50 6 454 282114] 

5 

Thus, the problem comprises dividing up the elements of Wmto seven groups, 
such that the sum of the elements in each group is close to but does not exceed 
811. We will now consider how, in the preferred embodiment, any fitness 
value and fitness vector is calculated for an individual or trial solution which is 
10 defined by an integer string b, of length 21, as follows: 

[4 53 7744 151267 5 623637 4] 

First, the twenty one integers are mapped to a trial phenotypic solution Mi 
15 matrix. Mi has seven columns, one for each of the seven possible groups, as 
follows: 



Group 


1 


2 


3 


4 


5 


6 


7 




310 


714 


170 


12 


12 


144 


176 


Mj 


497 


688 


50 


231 


380 


660 


216 




0 


0 


454 


266 


153 


6 


128 




0 


0 


0 


114 


0 


0 


282 


Sum 


807 


1402 


674 


623 


545 


810 


802 


Phenotype 












s 2 


S3 



schema 



20 



As may be evident, Mi is filled with the values taken from W, those values 
going into the columns defined by the integers of the string b. Thus, the first 



10 



WO 03/079288 PCT/GB03/00998 

w 9 

value (12) within Wgoes into column 4, the second value (12) goes into column 
5, the third value (170) goes into column 3, and so on. If a column already has 
a number in it, the next number is simply placed below. 

As shown above, the columns are then summed, and each sum is tested against 
the problem constraint conditions to see whether it complies or not. The 
problem constraint conditions will of course vary from problem to problem, but 
in this particular example the conditions are that the phenotype Mi for the string 
b should satisfy the two following conditions: 

sum(*j) < p = Sl\ 

Nc 

and £ (/> - sum («j) ) < p x m- sum (W) = 811 x 7 -5663 = 14 



forj=l,2,3 

where N c is the number of phenotype schemata (here 3) included in the 
1 5 phenotype My, and Sj is the y-th phenotypic schemata. 

In the table above, s u s 2 and s 3 are the phenotype schemata which satisfy the 
constraint conditions. Since there are three of those, the fitness value of the 
string b is 3 . The positions of the schemata mean that the fitness vector for b is 
20 [10 00 0 1 1]. 

More formally: 



• Assume that there are a total q of different phenotypic schema s u s 2> ... , s q 
25 searched in a population, the phenotype schemata vector M = [s u $2, ••• » S <H 
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is built up as a database of all searched phenotypic schemata in the entire 
population. 

• The phenotype of individual j can be mapped into its relative fitness binary 
vector /vQ as a binary string in relation to the schemata vector M according 
to the following principle: 

P ( j ) = [bi t>2 \] ; and bi ( z = 1, q) is a binary bit , 

(1 if the phenotype of individual j includes the schema Si; 
0 if the phenotype of individual j does not include the schema si 

To take a final example, the table below shows the phenotypic schema, fitness 
value and corresponding fitness vector for four individuals jj 9 j 2 , js> J*4 9 and 5 
searched phenotypic schemata, namely a, b, c, d, e: 



Individual 


Phenotypic 
schema 


Fitness value 


Fitness vector 


h 


fa,d > 


2 


ri. o, o, i,o) 




(b,c,d,e> 


4 


ro, i, l, i,n 


h 


{a,c,d,e} 


4 


( 1, 0, 1, 1, 1 ) 


h 


{a,b,c,d,e} 


5 


fi,i, i, i,i) 



Once the fitness value and fitness vector have been calculated for each of the 
individuals within the population, a selection process takes place to choose 
which of those individuals are to be retained and which discarded. The 
selection process is carried out in two stages: 



1. First, a reserved sub-population of the best individuals is kept. 
Individuals are chosen for this sub-population on the basis of their fitness 
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vectors (and possibly their fitness values as well). One copy only of 
each individual is retained, thereby ensuring that none of the best 
individuals within the previous generation can be lost. Typically, the 
size of this sub-population is limited to around 20% of the total 
5 population size. 

2. Next, the remaining 80% of the spaces are filled, using a so-called 
roulette wheel section method on all of the existing individuals. Using 
this method, offspring for the next generation are created on the basis of 
10 the fitness values of the individuals in the current generation. The 

higher the fitness value of an individual, the more offspring (copies) of 
that individual will be retained in the next generation. Preferably, the 
number of offspring is substantially proportional to the fitness value. 



15 In order to understand how the reserved 20% sub-population is filled, it is 
necessary to introduce the concept of Pareto Optimality: [Goldberg, D.E: 
Genetic Algorithms in Search Optimization, and Machine Learning, Reading, 
MA: Addison Wesley, 1989]. Using this concept, it may be stated that a vector 
x is partially less than j% symbolically x < p y , when the following conditions 

20 hold: 

jc<pj; <S> (V;)(x,. y / )A(3z)(x I <yO 



Under these circumstances we say that vector x dominates vector y. If a vector 
is not dominated by any other, we define it as nondominated. 

25 

Turning back to the previous table, giving fitness vectors for the individuals j h 
]2> js> j4» we can sa Y that individual j\ is dominated by both individual j 3 andy 4 , 
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but individual j 2 and j 3 are not dominated with each other. Individual z 4 is not 
dominated by any other individuals in the population: it is a nondominated 
individual in this population. 



5 Different non-dominated individuals should be kept in the next generation after 
selection. Actually, most of these non-dominated individuals contain few 
schemata; they are inferior individuals with relatively smaller absolute fitness 
values, and might not continue to evolve. So these non-dominated individuals 
are ranked according to their absolute fitness value: non-dominated individuals 

10 with greatest absolute fitness value Ma so far are ranked as 'first class' 
individuals, and others with absolute fitness value (Ma-k) (for positive integer 
k=l, 2,...; and k<Ma) are ranked as (Jc+1) class individuals. The value of £ is 
determined by the initial population size and the number of non-dominated 
individuals; the greater the value of k, the more genetic diversity will be 

15 preserved, and the more computation complexity will be required. In this case 
study, k is chosen as 1, which means mat only first-class non-dominated 
individuals searched so far are guaranteed to be preserved in the next generation 
in order to maintain suitable genetic diversity and accelerate convergence 
speed. 

20 

The number of non-dominated individuals is preferably fixed to a portion of 
population size, otherwise no space is available for new entries in the 
population. The size of non-dominated individuals is preferably set as 20% of 
population size. 

25 

Once the reserved sub-population has been determined, the remaining 80% is 
filled, as mentioned above, using roulette wheel selection: Goldberg, D.E.: 



* 
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Genetic Algorithms in Search, Optimization, and Machine Learning, Reading, 
MA: Addison Wesley, 1989. In this selection, bias optimum is used as the 
reproduction operator. The number of expected individuals is given by: 



integer ( FitnessValueof Individual x p opulation size ) 
& ^FitnessValueqflndividuals 



Once the new population has been chosen, crossover and mutation is applied to 
it. 



10 The basic forms of crossover are single-point, double-point, multi-point and 
uniform crossover. Single-point crossover creates two offspring from two 
parents. The parents are randomly selected from population, the crossover site, 
c x , is selected at random, and two offspring are made by both concatenating the 
bits that precede c x in the first parent with those follow c x in the second parent 

1 5 and performing the obverse operation. 

For multi-point crossover, multiple crossover points are chosen at random with 
no duplicates and sorted in ascending order. Then, the variables between 
successive crossover points are exchanged between the two parents to produce 
20 two new offspring. The section between the first variable and the first crossover 
point is kept between individuals. The following example illustrates this 
process. 

Consider the following two individuals with 1 1 binary variables each: 

25 

individual 1 01110011010 
individual 2 10 10 110 0 10 1 



♦ 
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cross points =3 3 6 10 
After crossover the new individuals are created: 

5 

offspring 1 0 1 1| 0 1 1| 1 1 0 1| 1 
offspring 2 1 0 1| 1 0 0| 0 0 1 0| 0 

The disruptive nature of multi-point crossover appears to encourage the 
10 exploration of the search space, rather than favouring the convergence to highly 
fit individuals early in the search, thus making the search more robust. 

In uniform crossover, a crossover mask which is a string of binary bits with the 
size of the chromosome, is generated randomly, the value of each bit in the 
15 mask determines for each corresponding bit in a child, which parent it will 
inherit that bit from. An example is given to illustrate the process. 



Parent 1: 2 5 1 4 3 6 

Parent 2: 4 6 2 1 5 3 

20 Mask: 0 10 0 11 

Child 1: 2 6 1 4 5 3 

Child 2: 4 5 2 1 3 6 



25 



Mutation is used as a secondary operator which is randomly applied with small 
probability. It is a safety policy to be used sparingly with reproduction and 
crossover against premature convergence to a local optimum. 



10 
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Two mutation operations are used here. One is "swap" where the two integers at 
two randomly selected distinct locations of the chromosome (individual) are 
exchanged. The other is "jump" where the integer in a randomly selected bit of 
the chromosome jumps to another feasible integer value. 

The mutation probability represents the frequency of applying the mutation 
operator. The mutation rate for jump method is determined by the following 
condition 

1 



m u — 



bitlength 



In addition to the previously mentioned techniques, further improvements can 
be achieved by including the following refinements: 

• Elitist selection strategy: the best individual is always kept in the next 
15 generation (see Eshelman. L.J: The CHC adaptive search algorithms: How 

to have safe search when engaging in non-traditional genetic 
recombination. Foundation of Genetic Algorithms. Gregory J.E. Rawlins 
(Eds), Morgan Kaufmann Publishers, pp.265-283, 1991). Also, the 
stochastic universal sampling method is used to reduce sampling bias (see 
20 Baker. J.E: Reducing bias and inefficiency in the selection algorithm. 

Proceedings of the Second International Conference on Genetic Algorithms, 
ppl4-19, 1987). 



25 



A highly disruptive crossover method (Eshelman, op cit) is used to increase 
the GA capability to exploring a broader search space to find prospective 
zones which may lead to a new fitness peak. 
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The method of the preferred embodiment may be implemented, without further 
inventive input, by a skilled programmer on the basis of the description given 
above. The algorithm may be conveniently implemented along the lines of the 
following pseudocode: 

Initialisation of population 
Generation = 0 

While termination condition not satisfied DO 

Construct phenotypic schemata vector M 9 fitness vector fv, through 

phenotypic cluster analysis 
Calculate the fitness value of all individuals 

Select suitable non-dominated individuals to form a sub-population whose 

size is a portion of overall population 
Select offspring to form new population using Roulette wheel selection 

method 

Apply highly uniform crossover and mutation to the new population 
Put all the sub-population into the new population 
Generation = Generation + 1 
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CLAIMS: 



1 . A method of combinatorial multimodal optimisation for finding multiple 
optimal ways of dividing a set W of n values into m groups, such that each of 
5 the groups satisfies a respective constraint condition, the method comprising: 



(a) defining an initial population of individuals, each representative 
of atrial solution; 

(b) calculating for each individual a fitness vector indicative of 
10 whether the constraint condition for each group has been satisfied; 

(c) selecting a plurality of individuals for the next generation in 
dependence upon their respective fitness vectors; 

(d) creating a new population including the selected individuals; and 

(e) repeating steps (b) to (d) until the population stabilizes, the 
15 individuals of the stable population representing multiple optional 

ways of dividing the set W. 



2. A method as claimed in claim 1 in which the fitness vector is of length 
m, each element in the fitness vector being indicative of whether the constraint 

20 condition of a corresponding one of the m groups has been satisfied. 

3. A method as claimed in claim 2 in which the fitness vector comprises m 
bits, each bit being indicative of whether the constraint condition of a 
corresponding one of the m groups has been satisfied. 

25 

4. A method as claimed in claim 1 including calculating a fitness value for 
each individual. 
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5. A method as claimed in claim 4 when dependent upon claim 3 in which 
the fitness value comprises the sum of the bits in the fitness vector. 

5 6. A method as claimed in any one of the preceding claims including 
reserving a proportion of the new population for individuals selected at step (c). 

7. A method as claimed 6 in which a non-reserved proportion of the new 
population is generated using a Roulette wheel selection method. 

10 

8. A method as claimed in any one of the preceding claims in which step 
(c) comprises selecting non-dominated individuals using the criteria of Pareto 
optimality. 

15 9. A method as claimed in claim 8 when dependent upon claim 4 including 
ranking non-dominated individuals by fitness value, and selecting from the 
ranked list. 

10. A method as claimed in claim 9 in which only non-dominated 
20 individuals with greatest fitness value may be selected at step (c). 

11. A method as claimed in claim 4 in which step (c) comprises selecting 
individuals in dependence upon both their respective fitness vectors and their 
respective fitness values. 
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12. A method as claimed in any one of the preceding claims in which 
crossover and mutation are applied at step (d) to at least some individuals in the 
new population. 

5 13. A method as claimed in any one of the preceding claims in which step 
(c) comprises selecting no more than one individual for each unique fitness 
vector. 

14. A method of distributing a plurality of tasks between a plurality of 
10 devices connected together to form a network, wherein each device has an 
associated constraint on the amount of tasks that it can perform per unit of time, 
the method comprising: 

(a) generating a plurality of trial solution allocations to form an initial 
15 population of allocations; 

(b) calculating for each allocation a fitness vector indicative of 
whether the constraint condition for each device has been satisfied; 

(c) selecting a plurality of allocations for inclusion in the next 
generation of allocations in dependence upon their respective fitness vectors; 

20 (d) creating the next generation of allocations by including the 

allocations selected in step (c) together with new allocations each of which is 
formed from a combination of two or more of the allocations selected in step 

(c); 

(e) repeating steps (b) to (d) until the population stabilizes; and 
25 (f) allocating the tasks among the devices according to one of the 

allocations included in the stabilized population. 



4 

r 



WO 03/079288 



22 



PCT/GB03/00998 



15. A method as claimed in claim 14 wherein the devices are processors 
within a multi-processor computer system. 

16. A method as claimed in either of claims 14 or 15 wherein the devices are 
5 computers within a computer network. 

17. A method as claimed in claim 14 wherein the devices are routers and the 
tasks are estimated volumes of traffic to be routed through the routers within a 
data network, and wherein the allocations are used to form a routing strategy. 

10 

18. A method as claimed in any one of claims 14 to 17 in which step (c) 
comprises selecting non-dominated allocations using the criteria of Pareto 
optimality of the associated fitness vectors. 

15 19. A method as claimed in any one of the preceding claims in which new 
allocations are formed in step (d) by performing crossover operations in respect 
of groups of two or more of the allocations selected in step (c). 

20. A method as claimed in any of claims 14 to 19 in which mutation 
20 operations are applied to one or more of the new allocations formed in step (d) 

according to a predetermined probability of each new allocation being mutated. 

21. A computer program for carrying out the steps of any one of the 
preceding claims. 

25 

22. A carrier medium carrying the computer program of claim 21 . 
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23. A system comprising a plurality of devices connected together to form a 
network, wherein each device has an associated constraint on the amount of 
tasks that it can perform per unit of time, the system including means for 
allocating a plurality of tasks among the devices, the allocation means 
5 comprising: 



(a) means for generating a plurality of trial solution allocations to 
form an initial population of allocations; 

(b) means for calculating for each allocation a fitness vector 
10 indicative of whether the constraint condition for each device has been satisfied; 

(c) means for selecting a plurality of allocations for inclusion in the 
next generation of allocations in dependence upon their respective fitness 
vectors; 

(d) means for creating the next generation of allocations by including 
15 the allocations selected in step (c) together with new allocations each of which 

is formed from a combination of two or more of the allocations selected in step 

(c); 

(e) means for repeating steps (b) to (d) until the population stabilizes; 

and 

20 (f) means for allocating the tasks among the devices according to one 

of the allocations included in the stabilized population. 



